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We analyze the full counting statistics of charge transfer through a quantum dot in the Kondo 
regime, when coupled to an arbitrary number of terminals N. At the unitary Kondo fixed point 
and for N > 2 we recover distinct anticorrelations of currents in concurring transport channels, 
which are related to the fermionic Hanbury Brown and Twiss (HBT) antibunching. This effect 
weakens as one moves away from the fixed point. Furthermore, we identify a special class of current 
correlations that are due entirely to the virtual polarization of the Kondo singlet. These can be 
used for extracting information on the parameters of the underlying Fermi-liquid model. 

PACS numbers: 72.10.Fk, 73.63.-b, 73.23.-b, 72.15.Qm 



Since its discovery the Kondo effect remains one of 
the most active fields in condensed matter physics 0. 
The last decade witnessed a very important achievement 
from the experimental point of view: the creation of a 
fully adjustable artificial Kondo impurity 0,0, Qt< This 
enables one, in principle, to measure almost any trans- 
port property of the system in a wide parameter range. 
One impressive example is the Kondo resonance splitting, 
which occurs under nonequilibrium conditions when the 
impurity is simultaneously coupled to two terminals held 
at different chemical potentials 0, 0, • It turned out 
that the best way to explore this effect is to introduce 
a third electrode and measure the tunnelling density of 
states 0, 0. This fact triggered a number of different 
studies of the Kondo effect in multiterminal geometries 

Apart from peak splitting measurements, artificial 
Kondo impurities are also interesting for another reason. 
As was pointed out in early works on transport in mul- 
titerminal systems, they possess some unique transport 
properties due to a possibility of having multiple con- 
current transport channels. In particular, the fermionic 
nature of charge carriers leads to antibunching of indi- 
vidual transport events, also called fermionic Hanbury 
Brown and Twiss (HBT) correlations |14L Hq . This re- 
veals itself as a negative sign of the current-current cor- 
relation function between concurring transport channels. 
For the simplest three-terminal geometry this effect has 
been confirmed experimentally in pH hj. Thus far 
little has been done to investigate these effects in the 
Kondo setup. The most advanced work is 0], where the 
current-current correlations have been analysed for the 
three-terminal geometry using the slave-boson approxi- 
mation. To the best of our knowledge, nothing is known 
about the correlations of higher orders or close to the 
fixed point. With this Letter we would like to close this 
gap by addressing the full counting statistics (FCS) of 
such setups near the unitary limit. 



We start with the following Hamiltonian, 
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where the first sum describes the dynamics of the elec- 
trons H [ipja] = J2 P crj( e p ~ V 3)% a j'>Pp<T3 in tnc noninter- 
acting leads (We set e = h = ks = 1 throughout. The 
conductance quantum then becomes Go = 2e 2 /h = 1/tt). 
Here tppo-j is the annihilation operator for an electron 
with momentum p, spin a in the j'th terminal held at 
chemical potential Vj. The second part represents the 
coupling of the local spin degrees of freedom in the leads 
to the impurity spin S: J is the respective coupling con- 
stant, r denotes the set of Pauli matrices, h — nsgB is 
the product of Bohr's magneton /is, of the gyromagnetic 
ratio g and the magnetic field B applied to the impu- 
rity. Via simple unitary transformation to a set of new 
fermions Cj a this Hamiltonian can be mapped onto the 
single channel Kondo model, where only one of the new 
fields, Cicr, remains coupled to the impurity spin: 
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In the low energy sector, well below the Kondo temper- 
ature Tk, this fermion degree of freedom is hybridised 
with the impurity spin into a singlet state. The true fixed 
point of the system is just an TV-leg resonant model. The 
effective fixed-point Hamiltonian is simply H for c\ a , 
supplemented by voltage terms. The best possible trans- 
parency is then realized by the following scattering ma- 
trix ipj(0 + ) = SjfcVfc(0~) (0 ± denote in- and out-going 
plane waves), with 18] 

Sjj =(N- 2)/N , s jk = -2/N for j^k. (3) 

The central object in FCS calculations is the cumu- 
lant generating function (CGF) X W = Eq e iQA P(Q) 
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[l9| . where A = (Ai, Ajv) are the measuring fields and 
P(Q) is the probability that charges Q = (Qi, ■■■,Qn) 
are transferred across the respective channels during the 
waiting time T. The cumulants (irreducible moments) 
are then found from the respective derivatives of x(^)- 
With the knowledge of the scattering matrix (no interac- 
tions) the CGF is given by the multiterminal version of 
the Levitov-Lesovik formula 20] : 



For V2 = V 3 the leading temperature dependence of this 
correlation is again linear, 



T 

lnxo(A) = — / dto In det 

TT 



i + /(w)(*ty-i) , (4) 



where s r , 
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and the diagonal matrix 
f{uj) =diag(ni,n 2 , . . . , n^r) contains the Fermi distribu- 
tion functions of individual terminals n.j. One then easily 
verifies that the fixed point CGF is given by: 
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In the particular case of zero temperature and three ter- 
minals such that Vi > V2 > V 3 the above formula be- 
comes 
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(Vi - V 2 ) In 1 + - (V^ 1 ^ + e t{Xl - x ^ - 2) 
1+4 ^i(A 1 -A 3 ) +e i(A 2 -A 3 )_2^ J . 
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Here one recognises that the current pairs Ii_>3 & I\^2 
and /i_>3 & ^2^3 are mutually correlated, the corre- 
sponding fluctuations of charge differences in the respec- 
tive channels being 



{S(Qi - Qi) S(Qi - Q 3 )) = 
(6(Qi-Qs)5(Q 2 -Q 3 )) 
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T(Vx-V 2 ), 
T{V 2 -V 3 ). (7) 



The negative sign is the signature of fermion antibunch- 
ing, or HBT type correlations. On the other hand, the 
pair ii_>2 & I 2 — »3 is not entangled in any way, as those 
transport channels are not concurring. This is true only 
at T = 0, of course, its leading temperature dependence 
being 



(S(Qi-Q 2 )S(Q2-Q 3 )) 
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for T <C |min(V^ — Vj)\, i 7^ j. In contrast to cross 
correlations, the autocorrelation functions are positive, 
e. g. one finds for the current noise between the channels 
2 and 3 
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One can show that these cumulants vanish whenever 
V 2 = V 3 and T = even away from resonance and for a 
lead- asymmetric junction. 

The departure from the Kondo fixed point is described 
by the Nozieres Fermi- liquid theory 21]. The effective 
Hamiltonian contains then in addition to the free part a 
scattering term 
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with a dimensionless amplitude a and local density of 
states at the Fermi energy of the leads v (assumed for 
simplicity to be equal for all terminals), and an interac- 
tion: 
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/ 2 T, 



K 



(cJ T ci T - n )(c\ l c n - n ) , (12) 



with an amplitude <j>. Here c\ a = J2 P c ipo- and no stands 
for the particle density in the zero field and accounts for 
the background ion charge. The constant a is responsible 
for the impurity specific heat whereas the sum a + <p 
governs the magnetic susceptibility of the system. In the 
actual Kondo model the Fermi-liquid relation a = <f> = 1 
holds |2ll. 

The standard way to calculate the FCS is to introduce 
time-dependent counting fields A (nonzero only du ring 
the waiting time T) coupled to charge currents 0, l22| . 
In the alternative procedure such terms can be gauged 
away by a canonical transformation psl ], which is easily 
generalized to the multi-terminal geometry. This trans- 
formation is performed at the expense of having A depen- 
dent exponentials embedded into the interaction parts of 
the Hamiltonian via flsc[ci] + -ffj n t[ c i] — * #a[ci]. The 
CGF is then calculated according to the prescription [24| , 



X (A) = xo(A)xi(A) 



= Xo(A) ( T c exp 



dtH x {t) 
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where C denotes the K eldy sh contour and Tq is the time 
ordering operator on it [23 |. The evaluation of XiW can 
now be accomplished via the standard linked cluster ex- 
pansion in the couplings a, <fi using the Green's functions 
of the A-rotated c\ a operator: 



2 N 
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where t± are the standard Pauli matrices, tq is the unity 
matrix and n JCT stand for t 
in the respective terminal. 



matrix and rij„ stand for the Fermi distribution function 
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To the lowest (second) order there are three different 
contributions: the ones proportional to a 2 , tfi 2 , and a<j>. 
The latter cross term is nonzero only in finite field. The 
a 2 part emerges as 



S a \nxi = -T7 
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Evaluation of the interaction 4> 2 correction results in 
two different contributions. The first one is field inde- 
pendent and for the case of N — 3 terminals (a general 
result for arbitrary N exists but is quite lengthy) is found 
to be 
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where we have defined the function 
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On the other hand, for the field-dependent correction we 
obtain 
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Finally, the a<f> cross term, which is also nonzero in a 
finite field only, is given by 
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The full FCS contains all transport characteristics of 
the system. Here we would like to concentrate on two 
rather interesting issues. The first question we want to 
address is to what extent are the HBT type correlations 
affected away from the fixed point. In order to discuss 
this we again use the convention V\ > V 2 = V3. Then the 
transport channels 1 — > 2, 1 — > 3 become concurrent. The 
only term, still containing both differences of counting 
fields Ai — A2,3, can be found in Eq. ljT5l) . The respective 
correction to J7J at T = is given by 

(§(Qi - Q2) 8{Qi - Q 3 )) = ^^r(Vi- V 2 f . (18) 



We see that the residual Kondo interactions lead to a 
weakening of the HBT signature. The remarkable fact is 
that this correction only depends on </> as well as its cu- 
bic voltage behaviour (in contrast to the linear one of the 
conventional HBT) offers an opportunity to extract this 
effective parameter from the experimental data. How- 
ever, the original HBT background might complicate the 
measurements. 

Luckily, this difficulty can be avoided by measurement 
of different correlation functions. As we have already 
mentioned above, at the fixed point there is no entangle- 
ment between nonconcurrent channels. This is different 
in the Kondo system. For instance, under the above con- 
ditions there is a residual correlation between the particle 
numbers in channels 1 and 3, which is absent at the true 
unitary limit. Using the result l|15(l we immediately ver- 
ify that 



(S(Qi-Q 2 )6{Q 2 -Q 3 )) = - 
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Comparing this equation with the temperature depen- 
dence of the primary contribution JSJ we see that the in- 
terplay of noncquilibrium and polarization of the Kondo 
singlet amounts to residual correlations at T = 0. Sim- 
ilar fate meets the autocorrelation (noise) of the I 2 — >3 
current. While it is zero in the unitary limit, Eq. for 
V 2 = V3, it acquires a finite value away from the fixed 
point, 



(s 2 (Q2 
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Since in this limit there is no background from the leading 
order, we expect that these types of correlation would be 
a good means to access the parameter in future 

experiments. Interestingly, the finite magnetic field does 
not affect these conclusions. 

As the system under consideration is a Fermi liquid 
it is natural to ask whether there is a possibility to use 
this fact in order to derive a universal CGF valid also 
beyond the leading corrections to the unitary limit. In- 
deed, in the case of two terminals Oguri succeeded in 
deriving an expression for the current for a more general 
Anderson impurity model for arbitrary on-site repulsion 
U for small voltages [26| . In |25| a related formula for a 
CGF has been proposed, valid around the unitary limit 
as well as for small U [13 • It turns out that a similar pro- 
gram can be performed in the present multiterminal case. 
Indeed, the cross-voltage derivatives, d 2 /dV\dV2, of the 
bare Green's functions vanish. Therefore the method of 
Ref.[2i| applies directly. For the relative currents [the 
actual currents in the respective channels being given by 
A = (l/JV)Ej in] we And 
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where Tjy = NT = Nirpoj 2 , 7 is the tunnelling ampli- 
tude coupling the Anderson impurity level to the termi- 
nals with the local tunnelling density of states po. The 
objects Xe,o are the even/odd susceptibilities which are 
correlations of impurity population probabilities with the 



same/different spin orientations, respectively. This result 
is valid for small voltages Vi -C Tjy and is exact in in- 
teraction as the U dependence is now contained in Xe,o- 
These are essentially equilibrium quantities and can be 
found with other techniques, e. g. by Bethe ansatz. In 
the unitary limit the identification of coefficients is very 
simple: Xe = T^a/Tx and Xo = Tff(f>/Tjc. After that 
substitution one then recovers the results found from the 
CGF lfT3|l . We expect that also in Eq. J5DJ) and similar 
formulas the parameter (f> can be substituted by the exact 
susceptibility Xo, even though we do not have a formal 
proof of this statement. 

To summarize, we report analytical results for the 
charge transfer statistics for the Kondo impurity coupled 
to an arbitrary number of terminals. In the Kondo limit, 
when bias voltage, magnetic field and temperature are 
all smaller than the Kondo temperature, we use a gener- 
alization of the Nozieres Fermi-liquid theory to derive a 
perturbative expansion for the FCS in leading irrelevant 
operators. At the true fixed point we recover with our 
method the Hanbury Brown and Twiss anticorrelations 
between concurrent transport channels. While these turn 
out to weaken away from the unitary limit, new correla- 
tions emerge due to the virtual polarization of the Kondo 
singlet. The most promising of these is the low-T noise in 
the 2 — 3 channel of the N = 3 junction at V\ > V% = V3. 
Measurements of the latter would enable one to directly 
access the effective parameters of the Fermi liquid model 
in experiments similar to those of [Tiil Il7| . 
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